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Water  transport  through  the  gas  diffusion  media  of  a  proton  exchange  membrane  fuel  cell  (PEMFC) 
was  investigated  with  a  focus  on  the  role  of  the  microporous  layer  (MPL)  coated  on  the  cathode  gas 
diffusion  layer  (GDL).  The  capillary  pressure  of  the  MPL  and  GDL,  which  plays  a  significant  role  in  water 
transport,  is  derived  as  a  function  of  liquid  saturation  using  a  pore  size  distribution  (PSD)  model.  PSD 
functions  are  derived  with  parameters  that  are  determined  by  fitting  to  the  measured  total  PSD  data. 
Computed  relations  between  capillary  pressure  and  liquid  saturation  for  a  GDL  and  a  double-layered 
GDL  (GDL  +  MPL)  show  good  agreement  with  the  experimental  data  and  proposed  empirical  functions. 
To  investigate  the  role  of  the  MPL,  the  relationship  between  the  water  withdrawal  pressure  and  liquid 
saturation  are  derived  for  a  double-layered  GDL.  Water  transport  rates  and  cell  voltages  were  obtained  for 
various  feed  gas  humidity  using  a  two-dimensional  cell  model,  and  are  compared  with  the  experimental 
results.  The  calculated  results  for  the  net  drag  with  application  of  the  capillary  pressure  derived  from  the 
PSD  model  show  good  agreement  with  the  experimental  values.  Furthermore,  the  results  show  that  the 
effect  of  the  MPL  on  the  cell  output  voltage  is  significant  in  the  range  of  high  humidity  operation. 
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1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC)  has  attracted 
a  great  deal  of  attention  in  the  last  decade  as  a  promising 
high-efficiency  low-emission  power  source  for  both  mobile  and 
stationary  applications.  However,  there  are  still  many  problematic 
issues  that  must  be  overcome  for  this  technology  to  be  efficient 
and  practical.  One  major  problem  is  water  management,  which 
requires  careful  attention  during  operation  of  a  state-of-the-art 
PEMFC.  At  low  humidity,  the  proton  exchange  membrane  and  elec¬ 
trode  assembly  (MEA)  lose  water,  which  leads  to  a  rapid  increase 
in  ohmic  resistance.  Conversely,  if  too  much  liquid  water  is  present 
in  the  cell,  then  the  pores  in  the  electrodes  are  filled  with  water 
and  the  passage  of  reactant  gases  becomes  obstructed.  Therefore, 
the  cell  operation  conditions  and  MEA  components  have  to  be  well 
matched  in  order  to  avoid  membrane  dehydration  and  cathode 
flooding.  In  addition  to  experimental  approaches,  water  manage¬ 
ment  models  are  useful  to  achieve  understanding  of  the  processes 
that  govern  water  transport  to  assist  optimization  of  the  fuel  cell 
operating  conditions  and  relevant  electrode  structures. 

Water  is  generated  in  a  PEMFC  at  the  cathode,  in  addition  to 
which,  water  is  transported  from  the  anode  to  the  cathode  by 
electro-osmotic  drag.  As  water  builds  up  at  the  cathode,  back  dif¬ 
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fusion  occurs.  Water  is  transported  from  the  cathode  back  to  the 
anode,  because  the  concentration  of  water,  i.e.  the  activity,  is  higher 
at  the  cathode  than  that  at  the  anode.  The  total  water  flux  occurs 
by  diffusion  through  the  gas  diffusion  media.  If  the  vapor  pressure 
reaches  saturation  pressure,  then  the  condensed  vapor  moves  as  a 
liquid.  In  such  a  case,  multiphase  flow  is  the  transport  mechanism 
for  water  flux.  Permeability  of  a  porous  medium  is  one  of  the  con¬ 
trolling  factors  that  affect  the  rate  at  which  fluids  travels  through 
pores.  Any  local  blockage  of  usually  open  pores  restricts  reactant 
flow  to  the  reaction  sites,  a  phenomenon  referred  to  as  flooding. 
Flow  of  the  liquid  water  in  a  porous  medium  is  driven  by  capillary 
action  generated  by  the  liquid  saturation  gradients. 

Capillary  pressure  increases  as  the  pore  radius  decreases;  there¬ 
fore,  it  is  important  to  determine  the  exact  pore  size  distribution 
of  the  porous  media  and  its  hydrophobic  or  hydrophilic  charac¬ 
ter.  These  parameters  are  usually  obtained  experimentally  and  are 
highly  dependent  on  the  material.  Various  measurements  of  cap¬ 
illary  pressure  have  been  reported  [1-6].  Kumbur  et  al.  presented 
a  modified  Leverett  function  to  estimate  the  capillary  pressure  as 
a  function  of  liquid  saturation  and  hydrophobic  additive  content 
based  on  experimental  data  [1-3].  Carbon  paper  (SGL24BC,  24CC, 
and  24DC)  and  carbon  cloth  (ELAT1200)  were  employed  as  gas 
diffusion  layers  (GDLs)  and  were  coated  on  one  side  with  a  micro¬ 
porous  layer  (MPL).  Fairweather  et  al.  conducted  experimental 
measurements  of  capillary  pressure  for  several  GDLs  (Toray  TGP- 
090  and  Avcarb  P75T)  after  hydrophobic  treatment  and  reported  a 
significant  hysteresis  between  the  liquid  intrusion  and  gas  intru- 
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Nomenclature 

List  of  symbols 

cij  activity  of  water  in  stream  j  (anode,  cathode) 

Cj  concentration  of  water  in  the  membrane,  mol  cm-3 

Dw  diffusion  coefficient  of  water,  cm2  s-1 
En  Nernst  loss,  V 

fr  k  fraction  of  pore  size  distribution,  k 
f  hi  fraction  of  hydrophilic  component 

F  Faraday  constant,  96,485  C  mol-1 

I  local  current  density,  A  cm-2 

i0j  exchange  current  density,  A  cm-2 

k  permeability  of  water,  m2 

krh  relative  permeability  of  water 

Mh2o  molar  weight  of  water 

nd  electro-osmotic  drag  coefficient 

Nh2o  molar  flux  of  water  (perpendicular  to  the  gas  diffu¬ 
sion  media),  mols-1  cm-2 
p  pressure,  kPa 

pc  capillary  pressure,  kPa 

R  gas  constant,  8.314  J  mol-1  K_1 

RH  relative  humidity 

r  pore  radius 

r0  k  characteristic  pore  radius  of  pore  size  distribution, 

k 

sk  spread  of  pore  size  distribution,  k 

si  liquid  saturation 

td  condensation  temperature,  °C 

T  temperature,  K 

V  cell  output  voltage,  V 

Voc  open  circuit  voltage,  V 

v(r)  normalized  volume  of  pore 

vp  pore  volume  of  porous  media 

x  coordinate  along  the  gas  channel 

y  coordinate  perpendicular  to  the  gas  diffusion  media 

Greek  symbols 

a  net  drag  of  water  per  proton 

/3k  coefficient  used  in  Eqs.  (6)  and  (7);  1  for  HI  pores 
and  -1  for  HO  pores 
y  surface  tension,  N  nrr1 

8  thickness,  cm 

s  porosity  of  the  gas  diffusion  media 

r\  overpotential,  V 

6C  water  contact  angle 

A  water  content  in  the  membrane 

r  tortuosity  factor 

§  stoichiometry 

a  conductivity,  cm-1 

v  viscosity,  Pa  s 

Subscripts 
a  anode 

c  cathode 

HI  hydrophilic 

HO  hydrophobic 

k  component  of  pore  size  distribution 

l  liquid 

nw  non-wetting 

sat  saturation 

total  total 

v  vapor 


sion  curves  [ 5 ].  Gostick  et  al.  also  reported  hysteresis  characteristics 
in  capillary  pressure  measurements  for  several  GDLs  (Toray  TGP- 
090  and  SGL10BA)  [6]. 

Weber  et  al.  have  developed  a  model  to  calculate  relevant 
parameters,  such  as  liquid  saturation,  permeability  and  the  aver¬ 
age  Knudsen  radius,  from  capillary  pressure  measurements  and 
structural  properties  such  as  the  pore-size  distribution  and  water 
contact  angle  [7]. 

There  have  been  several  reports  regarding  the  effect  of  a  MPL 
on  a  GDL  [7-1 1  ].  Weber  et  al.  insisted  that  a  MPL  enhances  fuel  cell 
performance,  due  to  an  increase  of  the  liquid  penetration  depth 
into  the  membrane  by  the  increased  pressure  [7,8].  Pasaogullari 
and  Wang  modeled  the  two-phase  (water  vapor  and  liquid  water) 
flow  and  transport  in  the  air  cathode  of  a  PEMFC  and  concluded 
that  capillary  action  is  a  dominant  transport  process  to  remove 
water  inside  the  two-phase  zone  [9].  Karan  et  al.  reported  that 
there  was  no  statistically  significant  difference  in  the  experimen¬ 
tal  net  drag  for  cells  with  MPLs  compared  those  without  MPLs 
[10,11]. 

In  this  work,  we  present  a  method  for  computing  the  capillary 
pressure  of  both  MPL  and  GDL  as  a  function  of  liquid  saturation 
using  the  pore  size  distribution  (PSD)  model.  The  model  derives  PSD 
functions  of  which  the  parameters  are  determined  by  fitting  to  the 
measured  total  PSD  data.  The  second  aim  is  to  investigate  the  role 
of  the  MPL  and  the  effect  of  capillary  pressure  on  water  transport 
and  cell  performance  based  on  the  PSD  model.  The  computation 
results  are  compared  with  experimental  water  transport  data. 

2.  Water  transport  characteristics  by  the  PSD  model 

2.1.  PSD  model 


Water  transport  mechanisms  in  gas  diffusion  media  and  catalyst 
layers  are  closely  dependent  on  the  pore  size  distribution  (PSD) 
of  the  media/layer.  Gas  diffusion  media  are  composite  structures 
that  exhibit  either  hydrophilic  or  hydrophobic  properties,  or  both. 
The  water  contact  angle  (0C)  for  composite  media  is  expressed  as 
follows: 


0C  =  cos  1(Jhi  *  cos  Ohi  +  (1  -fm)  *  cos  Oho), 


(1) 


where /h/  is  the  fraction  of  hydrophobic  characteristic,  and  0Hi  and 
0Ho  are  the  water  contact  angles  of  the  hydrophilic  and  hydropho¬ 
bic  components,  respectively.  Scanning  electron  microscopy  (SEM) 
images  of  gas  diffusion  media  show  that  pores  coated  with  poly- 
tetrafluoroethylene  (PTFE)  consist  of  both  a  larger  hydrophobic 
portion  and  a  smaller  hydrophilic  portion.  It  has  been  reported 
that  the  fraction  of  hydrophobic  pores  in  bare  carbon  paper  ranges 
between  2%  (SGL10BA)  and  15%  (TGPH  120)  [12].  Therefore,  it 
is  assumed  that  the  water  contact  angle  of  a  hydrophobic  pore 
is  approximated  by  the  combination  of  a  90%  hydrophobic  com¬ 
ponent  and  a  10%  hydrophilic  component  using  Eq.  (1).  On  the 
contrary,  a  hydrophilic  pore  is  approximated  by  the  combination 
of  a  10%  hydrophobic  component  and  a  90%  hydrophilic  compo¬ 
nent.  We  have  assumed  contact  angles  for  bare  carbon  fiber  and 
PTFE  of  80°  and  110°,  respectively  [7,13].  Surfacial  contact  angles 
of  hydrophobic  and  hydrophilic  pores  were  calculated  using  Eq.  (1 ) 
as  83°  and  107°,  respectively,  which  were  in  agreement  with  the 
measured  contact  angles  for  bare  and  hydrophobic-treated  carbon 
papers  [13]. 

The  PSD  of  the  gas  diffusion  media  is  expressed  by  the  summa¬ 
tion  of  log  normal  distributions  originally  described  by  Weber  et  al. 
[7], 


v(r)  = 

k 


1 

rskV2jr 


exp 


1  /  In  r  -  In  r0,fc 

2V  skV2 


2 


(2) 
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Table  1 

Integration  limits  of  the  PSD  model  and  coefficients  in  Eqs.  (6)  and  (7). 
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Pores  filled 
with  water 

Integration 

limits 

Coefficient  in 

Eqs.  (6)  and  (7) 

Hydrophilic  (HI) 
Pc<  0 

r<rCiHI 

0  to  rCtHi 

Pc>  0 

All  pores 

0  to  oo 

Hydrophobic  (HO) 
Pc<  0 

_ 

_ 

Pc>  0 

rc,HO  <  r 

rc,HO  <  oo 

where  v{r)  is  the  normalized  volume  of  pore  radius  r,  rok ,  and  sk  are 
the  characteristic  pore  radius  and  spread  of  distribution  /<,  respec¬ 
tively.  fr  k  is  the  fraction  of  distribution  k  in  the  total  distribution. 
These  values  are  determined  from  fitting  with  the  experimental 
PSD  data.  A  pore  of  given  radius  r  is  determined  as  being  either  filled 
or  empty  at  a  capillary  pressure  pc  according  to  the  corresponding 
pore  radius  rc  defined  by  Eq.  (3), 

fc  =  -—  COS  —  .  (3) 

Pc 


Liquid  saturation  si  and  relative  liquid  permeability  krh,  are 
important  factors  used  to  describe  the  two-phase  transport  in  the 
gas  diffusion  media,  si  is  defined  as  the  volume  fraction  of  the  vol¬ 
ume  occupied  by  the  liquid  phase  vpi,  divided  by  the  total  void 
volume  of  the  porous  media  vp<totai, 


Vp,  total 


Jrrcfv(r)dr 
Jfv(r)dr  ' 


(4) 


krh  is  defined  as  the  ratio  of  the  effective  liquid  permeability 
through  the  pores  occupied  by  the  liquid  phase  kh  to  the  perme¬ 
ability  at  complete  saturation  ksat • 


ffj  r2v(r)dr 
ffr2v(r)dr 


(5) 


sl  and  krh  can  be  calculated  by  integrating  and  summing  for  each 
region  where  liquid  water  is  present  using  Eqs.  (2),  (4)  and  (5).  The 
integration  regions  (rc  i  to  rc  2)  are  different,  depending  on  whether 
a  pore  is  hydrophilic  or  hydrophobic  and  are  shown  in  Table  1.  The 
corresponding  pore  radii  rc  at  pc  are  calculated  from  Eq.  (3)  and 
defined  as  rc<HI  and  rc<H0  for  hydrophilic  and  hydrophobic  pores, 
respectively.  Eqs.  (6)  and  (7)  are  then  derived: 


s'=z¥{i+^( 

krh  =  J2f-f\l+Pkerf 


In  rc  h  -  In  r0 


SkV2 

In  rc  h  -  In  r0Jf 
skV2 


*)}• 


-SkV2)  , 


(6) 

(7) 


where  erf  represents  the  error  function.  The  coefficient  fik  comes 
from  the  difference  of  integration  limits  and  is  defined  as  1  for 
hydrophilic  pores  and  -1  for  hydrophobic  pores. 


2.2.  Relation  between  the  capillary  pressure  and  liquid  saturation 

Several  measured  PSD  curves  have  been  reported  for  commer¬ 
cial  carbon  paper  and  cloth  after  hydrophobic  treatment,  usually 
with  MPLs  [1-3,12,14].  Kumbur  et  al.  [1]  measured  hydrophilic 
and  total  PSD  curves  for  carbon  paper  and  cloth  using  the  standard 
porosimetry  technique  developed  by  Porotech.  We  have  applied 
the  PSD  function  of  Eq.  (2)  to  the  PSD  curve  measured  for  a  dou¬ 
ble  layered  gas  diffusion  medium  reported  by  Kumbur  et  al.  (Fig.  3 
in  [1]).  The  total  and  hydrophilic  PSD  curves  are  well  represented 
by  two  hydrophilic  PSDs  (HI1  and  HI2)  and  four  hydrophobic  PSDs 


Fig.  1.  Computed  pore  distribution  for  gas  diffusion  media.  ♦:  total  pores;  ■: 
hydrophilic  pores. 


(H01-H04),  as  shown  in  Fig.  1.  The  parameters  of  these  PSDs  are 
summarized  in  Table  2. 

The  relationship  between  capillary  pressure  and  liquid  satu¬ 
ration  can  be  derived  by  applying  Eqs.  (3)  and  (6)  to  the  PSD 
functions.  The  derivation  process  is  as  follows:  firstly,  a  pressure 
is  assumed  and  the  liquid  volume  filled  in  six  pore  distributions  are 
calculated  and  liquid  saturation  is  derived  by  summation.  The  cal¬ 
culation  is  repeated  assuming  different  pressures,  and  the  results 
are  shown  in  Fig.  2.  Negative  and  positive  capillary  pressures  cor¬ 
respond  to  hydrophilic  and  hydrophobic  pores,  respectively,  and 
capillary  pressure  is  zero  at  sl  =  0.222.  The  capillary  pressure  varies 
stepwise  in  the  curve,  which  reflects  that  the  total  PSD  is  composed 
of  discrete  PSDs.  The  capillary  pressure  for  an  HOI  micropore  is  out 
of  the  range  of  this  figure,  because  the  capillary  pressure  range  is 
more  than  100  kPa.  This  curve  can  be  expressed  approximately  by 
a  multi-power  equation,  which  is  represented  as  the  solid  line  in 
Fig.  2: 

Pc  (kPa)  =  -75.723s/5  +  440.53s/4  -  495.54s/3  +  161.14s/2 

+  24.358s/ -7.8799.  (8a) 


Kumbur  et  al.  proposed  modified  Leverett  functions  for  a  GDL 
with  MPL  to  estimate  the  capillary  pressure  as  a  function  of  non¬ 
wetting  liquid  saturation  and  hydrophobic  additive  content  for  the 
hydrophobic  pores  [1  ].  The  functions  were  approximated  with  the 
experimental  capillary  pressure  data  for  carbon  paper  with  PTFE 
contents  ranging  from  0  to  20  wt%  and  with  the  MPL.  The  calculated 
Eq.  (8a)  is  comparable  with  these  empirical  functions.  The  modi- 


Table  2 

Structural  properties  of  the  gas  diffusion  media. 


GDL 

Thickness  d  0.019 

Hydrophobic  contact  angle  Oho  107 

Hydrophilic  contact  angle  0Hi  83 

Porosity  e  0.8 

Tortuosity  factor  t  5 

MPL 

Thickness  8  0.005 

Hydrophobic  contact  angle  0Ho  107 

Porosity  e  0.5 


Pore  size  distribution  characteristics 


Pore 

Characteristic  radius 

Spread 

Fraction 

k 

1~ok 

sk 

frk 

HOI 

0.05 

0.25 

0.21 

H02 

1 

0.5 

0.04 

H03 

4 

0.15 

0.36 

H04 

9 

0.1 

0.18 

HI1 

3 

0.1 

0.14 

HI2 

15 

0.2 

0.07 
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Liquid  saturation 

Fig.  2.  Relation  between  capillary  pressure  and  liquid  saturation  for  gas  diffusion 
media.  ♦:  estimation  for  a  double-layered  GDL,  solid  line  is  the  approximation  from 
Eq.  (8a);  ■:  estimation  for  a  GDL,  dotted  line  is  the  approximation  from  Eq.  (8b). 

tied  Leverett  functions  consist  of  three  different  expressions  that 
correspond  to  the  different  liquid  saturation  regions:  0  <  slnw  <  0.5, 
0.5  <  slnw  <  0.65,  and  0.65  <  slnw  <  1 .  The  slnw  values  correspond  to 
the  hydrophobic  region,  so  that  it  is  necessary  to  transform  slnw  =  0 
in  the  expression  to  si  =  0.222  from  our  calculation  for  comparison. 
The  modified  Leverett  function  (PTFE  5  wt%)  is  in  good  agreement 
with  the  results  calculated  from  Eq.  (8a)  for  sl<  0.59,  as  shown  in 
Fig.  3a. 

2.3.  Pressure  for  water  transport  in  gas  diffusion  media 

2.3. 1.  Gas  diffusion  media  ( GDL  and  MPL) 

A  highly  hydrophobic  layer  on  a  GDL  is  referred  to  as  a  MPL. 
A  hydrophobic  HOI  PSD  (characteristics  shown  in  Table  1)  con¬ 
tributes  the  main  part  of  the  PSD  for  the  MPL,  as  described  by 


0.0  0.2  0.4  0.6  0.8  1.0 

Liquid  saturation 

Fig.  3.  (a)  Comparison  of  the  simulation  results  for  capillary  pressure  as  a  function  of 
liquid  saturation  using  the  modified  Leverett  function  from  Ref.  [  1  ].  The  dotted  line  is 
the  approximation  from  Eq.  (8a)  and  the  solid  line  is  the  modified  Leverett  function, 
(b)  Comparison  of  the  computed  capillary  pressure  curves  with  experimental  data. 
♦,  0:  Fairweather  et  al.  (♦:  liquid  intrusion,  0:  gas  intrusion)  [5];  ■,□:  Gostick  et  al. 
(■:  liquid  intrusion,  □:  gas  intrusion)  [6].  The  dotted  line  is  the  approximation  from 
Eq.  (8b). 


Kumbur  et  al.  [1].  However,  it  is  not  clear  whether  the  other 
hydrophobic  PSD  (H02)  belongs  to  either  the  GDL  or  MPL.  SEM 
images  indicate  that  the  MPL  is  microporous  with  macro-cracks 
[13],  which  suggest  that  the  MPL  encompasses  both  hydrophobic 
HOI  and  H02  PSDs. 

The  validity  of  the  computed  PSD  characteristics  for  double¬ 
layered  gas  diffusion  media  (GDL  +  MPL)  can  be  evaluated  in  terms 
of  a  comparison  with  the  gas  permeability.  Assuming  that  Hagen- 
Poiseuille  flow  is  applicable  for  porous  media,  then  the  relationship 
between  permeability  I<  and  characteristic  pore  radius  r,  and  the 
porosity  s  is  determined  as  [15]: 


where  kk  is  the  Kozeny  constant,  which  represents  the  shape  factor 
multiplied  by  the  tortuosity  factor.  When  a  porous  medium  consists 
of  several  specific  pores,  then  each  permeability  should  be  evalu¬ 
ated  using  Eq.  (9)  and  summated.  The  following  permeabilities  are 
calculated  for  double-layered  and  bare  GDLs,  respectively,  assum¬ 
ing  that  the  pores  are  circular  (kk  =  2  for  circular  pores)  and  applying 
the  characteristic  pore  radius,  fraction  and  porosity  shown  in 
Table  1 :  I<dGDL  =  0.121  x  1 0~u  m2  and  I<GDL  =  3.74  x  1 0~u  m2.  These 
values  approximately  coincide  with  the  reported  permeabilities  for 
double  layered  and  bare  GDLs,  I<dGDL  =  0.07-0.58  x  10-12  m2  and 
I<gdl  =  8- 69-31  x  10-12m2,  respectively  [14],  which  also  suggests 
that  the  PSDs  of  the  MPL  consist  of  two  types  of  PSD:  HOI  and  H02. 

The  relationships  between  capillary  pressure  and  liquid  satu¬ 
ration  for  a  GDL  can  be  derived  by  applying  Eqs.  (3)  and  (6)  to 
the  PSD  functions  for  HI1,  HI2,  H03,  and  H04  using  the  same  pro¬ 
cess  described  in  Section  2.2.  These  relationships  can  be  expressed 
approximately  by  a  multi-power  equation  as  represented  by  the 
dotted  line  in  Fig.  2: 

Pc  (kPa)  =  -204.37sf5  +  604.99s/4  -  625.82 si3  +  249 si2 

-4.4181s/ -7.1964.  (8b) 

Several  experimental  data  for  pc  vs.  si  reported  by  Fairweather  et 
al.  (Toray  TGPH  090)  [5]  and  Gostick  et  al.  (SGL 10BA  with  5%  PTFE) 
[6]  are  also  plotted  in  Fig.  3b  for  comparison  with  the  calculated 
data.  The  capillary  pressure  calculated  using  Eq.  (8b)  shows  nega¬ 
tive  values  for  si  <  0.2  and  increases  gradually  up  to  approximately 
lOkPa  at  sl^  0.8.  These  characteristics  are  in  agreement  with  the 
measured  capillary  pressure  under  conditions  of  water  intrusion. 

2.3.2.  Relationship  between  the  pressure  for  water  transport  and 
liquid  saturation 

The  relationship  between  the  pressure  for  water  transport  and 
liquid  saturation  for  the  GDL  is  derived  using  PSD  functions  for  HI1 , 
HI2,  H03  and  H04,  and  that  for  the  MPL  is  derived  using  PSD  func¬ 
tions  for  HOI  and  H02,  and  the  results  are  shown  in  Fig.  4a  and  b, 
respectively.  In  Fig.  4a,  the  pressure  data  for  hydrophilic  pores  are 
plotted  as  positive  capillary  pressures,  because  the  necessary  pres¬ 
sure  for  the  withdrawal  of  liquid  water  from  the  porous  media  is 
positive.  The  relationships  for  a  GDL  are  expressed  approximately 
by  the  multi-power  Eqs.  (10a)  and  (10b),  where  s/0  is  the  liquid 
saturation  corresponding  to  zero  capillary  pressure.  Pores  corre¬ 
sponding  to  less  than  s/0  are  hydrophilic  and  are  always  filled  with 
water  under  steady  state  operation. 

In  Fig.  4b,  two  curves  are  shown  for  the  relationship  between 
pressure  for  water  withdrawal  and  liquid  saturation.  The  curves 
pc-s/( a)  and  pc-s/(b)  correspond  to  MPLs  that  consist  of  both 
hydrophobic  HOI  and  H02  PSD,  and  only  hydrophobic  HOI  PSD, 
respectively.  The  capillary  pressure  ofpc-s/(a)  is  much  lower  than 
that  ofpc-s/(b)  in  the  low  liquid  saturation  region,  due  to  the  pres¬ 
ence  of  macro-cracks  (H02).  Hereafter,  we  apply  the  pc-s/(a)  curve 
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Fig.  4.  (a)  Pressure  for  water  withdrawal  from  a  GDL  as  a  function  of  liquid  satu¬ 
ration.  ■  :  calculated  pressure  for  water  withdrawal.  The  dotted  line  represents  the 
approximation  from  Eq.  (10a),  and  the  solid  line  represents  the  approximation  from 
Eq.  (10b).  (b)  Pressure  for  water  transport  through  the  MPL  as  a  function  of  liquid 
saturation.  ■:  pc-sl( a),  ♦:  pc-s/(b).  The  solid  line  is  the  approximation  from  Eq.  (10c). 


as  representative  of  a  MPL.  The  relation  for  a  MPL  is  approximated 
by  Eq.  (10c): 

GDL: 

PUH0  =  1 182.6s/5  -  3585.4s(4+ 4163.2sl3  -  2311. Is/2  +627.62s( 
-63.873,  si  >  s/0,  (10a) 


P/w=-6652.2s/3+4950.3s/2  — 1232.4S/+103.58,  sl<sl0,  (10b) 
MPL: 

P,=  -  29,  269s/5+77,  716s(4-70,  085s/3+23,  329s(2-558.14s(.(10c) 

Almost  the  same  relation  is  derived  for  the  relative  permeabil¬ 
ity  constant  krh  and  the  capillary  pressure  P/,  i.e.,  the  relationship 
krh  =  si  exists. 

2.4.  Water  transport  by  capillary  pressure  in  gas  diffusion  media 


by  integrating  this  equation  in  all  regions  of  the  porous  media: 

p8t  7»S/2  .p 

nh2ojmh2o  /  dy=  krhS^dsl 
JO  J  sf 

=  [  °  krh(sl-sl0)2^fdsl 

Jsh  dsl 

+  [  2  krho(sl-sl0f^f^dsl,  (12) 

Jsi0  dsl 

where  8  is  the  thickness  of  the  GDL  or  MPL,  and  sf  and  sl2  are  the 
lower  and  upper  limits  of  liquid  saturation  for  liquid  flow.  Liquid 
saturation  sf  belongs  to  the  hydrophilic  region  and  sl2  belongs  to 
the  hydrophobic  region,  and  1  -  sl2  is  the  gas  diffusion  space.  There 
following  relationship  exists  between  sf  and  sl2, 

Pi,hi(sI  i )  =  Pi,Ho(sh)-  (13) 

P/  H/  and  P/  ho  are  multi-power  functions  of  si,  so  that  by  apply¬ 
ing  the  approximate  expression  krh  =  si.  Eqs.  (12)  and  (13)  are 
expressed  as  multi-power  functions  of  sf  and  sl2.  These  functions 
are  solved  by  Newton  method  to  obtain  sf  and  sl2. 

The  presence  of  liquid  water  in  gas  diffusion  media  obstructs 
the  diffusion  of  reactant  gases.  The  diffusion  space  for  cathode 
gas  is  reduced  from  1  to  the  fraction  1  -sl2.  The  partial  pressure 
drop  of  oxygen  in  the  GDL  becomes  larger  as  the  gas  diffusion 
space  is  reduced,  which  causes  an  increase  in  the  concentration 
polarization,  si  at  breakthrough  indicates  the  threshold  where  the 
movement  of  accumulated  liquid  water  through  the  GDL  is  initiated 
by  a  capillary  pressure-controlled  percolation  process.  The  accu¬ 
mulation  is  the  result  of  several  effects,  including  the  existence  of 
the  hydrophilic  pores  and  the  so-called  ink-bottle  effect.  In  the  ink- 
bottle  effect,  as  liquid  enters  a  widening  pore,  the  capillary  forces 
change  with  the  increasing  diameter,  which  can  trap  a  droplet  [16]. 
Therefore,  the  liquid  saturation  at  breakthrough  corresponds  to  the 
diffusion  space  occupied  by  liquid  water.  Gostick  et  al.  reported  liq¬ 
uid  saturation  at  breakthrough  ranging  between  0.14  and  0.20  for 
Toray  120C  and  060C  [17].  The  value  of  0.286  for  s/0  in  our  calcu¬ 
lation  (Fig.  4a)  is  almost  the  same  as  the  reported  liquid  saturation 
at  breakthrough. 

Calculation  of  liquid  saturation  gives  sf  =0.181  and  sl2  =  0.344 
for  the  GDL  at  the  baseline  operating  condition  shown  in  Table  2. 
The  MPL  consists  of  only  hydrophobic  pores;  therefore,  the  relation 
sl0  =  0  always  exists.  It  is  then  possible  to  obtain  sl2  by  solving  only 
Eq.  (12).  A  similar  calculation  for  the  MPL  leads  to  sl2  =  0.056.  The 
change  of  the  diffusion  space  as  a  function  of  the  water  flux  was  cal¬ 
culated  using  Eqs.  (10a)-(10c)  and  (12),  and  is  shown  in  Fig.  5.  The 
relationship  between  si  and  NH2q  /  shown  in  Fig.  5  is  approximated 
by  Eq.  (14): 

1  -s/2  =  0.7143  -231. 25N°26(™8.  (14) 


The  transport  of  liquid  water  from  a  cathode  catalyst  layer  to 
the  cathode  flow  channel  is  expressed  as  follows: 

m  n  n  krhS^K  .... 

Nh2o,/Mh2o  =  — - — VP/,  (11) 

where  the  relation  NH2o,/  =  /(I  +  2a)/ 2 F  exists  for  the  saturation 
condition  and  the  relation  5^  =  {si  -  s/0)2  is  the  correction  factor, 
which  reflects  the  separation  and  joining  of  liquid  flow  in  gas  dif¬ 
fusion  media,  as  derived  by  Weber  et  al.  [7].  The  hypothesis  is  that 
transfer  between  liquid  flow  paths  only  occurs  when  the  paths  are 
of  the  same  diameter.  The  right  side  term  in  Eq.  (11)  is  the  water 
flux  of  Darcy  flow  in  the  porous  media,  where  krhS^K  and  vP/  are 
the  relative  permeability  of  water  and  the  liquid  pressure  gradient 
in  the  GDL  or  MPL,  respectively.  The  following  relation  is  derived 


Fig.  5.  Reduction  of  cathode  reactant  gas  diffusion  space  as  a  function  of  water  flux. 
The  solid  line  is  the  approximation  from  Eq.  (12). 
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The  maximum  water  flux  (6.3  x  10-6molcm-2  s-1  in  Fig.  5)  cor¬ 
responds  to  a  current  density  as  large  as  1.0  A  cm-2.  The  diffusion 
space  decreases  from  0.714  to  0.550  for  this  water  flux  range.  The 
calculated  partial  pressure  drop  of  oxygen  in  the  GDL  using  the 
two-dimensional  cell  model  (model  formula  are  described  in  Sec¬ 
tion  3.1 )  is  approximately  1%  and  the  diffusion  space  is  0.645  under 
the  baseline  cell  operation  conditions.  Therefore,  a  diffusion  space 
of  0.550  at  a  current  density  of  1.0  A  cm-2  suggests  there  is  no 
significant  increase  in  the  concentration  polarization. 

3.  The  role  of  MPL  in  water  transport  from  the  anode  to 
cathode 

3.1.  Cell  model  for  simulation 

The  catalyst  layer  in  the  PEMFC  is  located  between  the  mem¬ 
brane  and  gas  diffusion  media,  and  water  is  generated  in  this 
layer.  There  may  be  a  pressure  change  in  the  catalyst  layer  when 
condensed  liquid  water  moves  to  the  membrane  or  cathode  gas  dif¬ 
fusion  media.  However,  the  pressure  change  in  the  catalyst  layer 
is  not  significant  compared  with  that  of  the  MPL  for  the  follow¬ 
ing  reasons:  (1)  the  short  distances  of  water  transport,  and  (2) 
larger  hydraulic  permeability  [8].  Therefore,  the  pressure  at  the 
interface  of  catalyst  layer/membrane  is  assumed  to  be  the  same 
as  that  at  the  catalyst  layer/MPL  interface.  The  pressure  of  the 
water  at  the  cathode  catalyst  layer/membrane  interface  is  higher 
than  the  saturation  pressure,  due  to  the  presence  of  the  hydropho¬ 
bic  MPL.  Therefore,  the  activity  of  water  at  the  cathode  catalyst 
layer/membrane  interface  is  increased  when  the  vapor  pressure 
reaches  the  saturation  pressure.  The  pressure  of  water  at  the  cath¬ 
ode  catalyst  layer/membrane  interface  is  expressed  by  Eq.  ( 1 5)  with 
addition  of  the  MPL  capillary  pressure  pc,mpl' 

P  =  Psat  +  Pc, MPL-  (15) 

By  applying  Eq.  (1 5)  to  the  water  activity  at  the  cathode  catalyst 
layer/membrane  interface  in  a  two-dimensional  cell  model  [18,19], 
it  is  possible  to  calculate  the  net  water  drag  a.  The  gas  diffusion  per¬ 
pendicular  to  the  MEA  (y-direction)  in  this  cell  model  is  expressed 
by  the  Stefan-Maxwell  equation  under  the  condition  of  no  nitro¬ 
gen  flux.  The  following  relationships  exist  between  the  molar  flux 
of  oxygen,  water  and  current  density  in  the  cathode  along  the  flow 
channel  (x-direction).  In  the  model,  the  cell  is  divided  into  eight 
sections  in  the  x  direction. 


£ 

II 

1 

(16) 

Nh2o,c  =  ^(1+2«). 

(17) 

Water  is  transported  in  the  membrane  by  electro-osmotic  drag 
caused  by  proton  transport  and  back  diffusion  due  to  the  activity 
(namely,  concentration)  gradient  of  water. 

N  „  IM  D  3 C(X)  J(x) 

Nh20 ,a  =  nd  F  Dw  dy  =a  F  . 

(18) 

The  electro-osmotic  drag  coefficient  n diffusion  coefficient  Dw, 
and  the  water  content  in  the  membrane  A,  are  shown  in  Eqs.  (19), 
(21)  and  (22),  respectively. 

nd =0.0049  +  2. 024a, -4.53a2 +4.09a3=l  .59  +  0.159(a,-l ), 

(19) 

Psat 

(20) 

Dw  =  1.1  x  10-6nd  exp  (2416  (^  -  1))  , 

(21) 

X  =  0.043  +  17.81a,  -  39.85a2  +  36a,3. 

(22) 

The  expressions  given  by  Springer  et  al.  for  nd,  and  X  are  applied 
[20].  The  expression  for  Dw  is  modified  (the  constant  is  twice  as 
large  as  that  of  the  previously  reported  value  [19]),  considering  a 
recently  reported  Fickian  diffusion  coefficient  that  was  experimen¬ 
tally  obtained  [21]. 

Eq.  (23)  shows  the  relationship  between  the  cell  output  voltage 
V,  cell  overpotential  p,  and  the  ionic  conductivity  of  the  membrane 
a. 


V  =  V0C- 


I(x)8m 

O 


{Pc  +  Pa)  ~  En, 


(23) 


where  Voc  and  En  are  the  open  circuit  voltage  and  Nernst  loss, 
respectively.  The  cathode  overpotential  pc ,  and  anode  overpoten¬ 
tial  pa,  are  given  using  the  exchange  current  densities  ioc  and  ioa , 
respectively,  as  follows: 


,  *RT  i 

Pc  =  1.5—  In 


(24a) 


Pa  =  0.5 


RT  /(x) 
F  ioaPu 2 


(24b) 


Based  on  the  discussion  in  Section  2.4,  another  modification  of 
the  cell  model  has  been  made  to  consider  the  reduction  of  the  gas 
diffusion  space,  expressed  by  Eq.  (14).  NH20  /(x)  is  calculated  by 
determining  the  saturated  water  vapor  from  the  water  flux  Nh2oM 
at  each  position  along  the  flow  channel. 


3.2.  Comparison  of  the  model  calculation  and  experimental  net 
drag 


3.2.1.  Experimental 

In  this  work,  we  present  a  method  for  computing  the  capillary 
pressure  characteristics  of  a  double  layered  GDL  using  the  reported 
PSD  data.  We  also  investigate  the  role  of  the  MPL  and  the  effect 
of  capillary  pressure  on  the  net  drag  using  the  derived  capillary 
pressure  characteristics  and  the  calculated  results  are  compared 
with  that  of  the  net  drag  measurement.  A  25  cm2  MEA  was  made 
of  a  Nation  112  membrane,  gas  diffusion  media  of  Toray  carbon 
paper  (TGH-090)  treated  with  10  wt%  PTFE  and  a  MPL  which  was 
made  of  carbon  black  and  PTFE  and  was  coated  on  one  side  of  the 
GDL  with  the  thickness  of  30  \xm.  Net  drag  was  derived  through  the 
measurement  of  the  water  condensed  at  the  inlet  and  the  exit  of 
the  test  cell.  Specification  of  the  test  cell  and  standard  experimental 
condition  are  summarized  in  Table  3.  Further  details  of  the  net  drag 
measurement  are  described  in  the  previous  report  [18]. 

3.2.2.  Comparison  of  the  calculation  and  experimental  on  the  net 
drag 

Although  it  is  well  known  that  properties  of  GDL  depend  on  its 
venders,  we  apply  the  results  on  the  capillary  characteristics  of  a 
double  layered  GDL(SGL24DC)  to  the  comparison  of  water  trans¬ 
port  characteristics  from  the  following  reasons:  (1 )  Although  there 


Table  3 

Specification  of  a  test  cell  and  standard  test  condition. 


Membrane  equivalent  weight 

1200 

Membrane  thickness  (cm) 

50  x  10-4 

Membrane  dry  density 

1.84 

Catalyst  layer  thickness  (cm) 

10 x 10“4 

GDL  thickness  (cm) 

300  x  10“4 

Current  density  (A cm-2) 

0.3 

Cell  temperature  (°C) 

80/70 

Fuel 

h2 

Cathode  reactant  gas 

Air 

Stoichiometry  of  anode  gas,  xa 

1.43 

Stoichiometry  of  cathode  gas,  xc 

2 

Pressure,  kPa 

101.3 

Humidity  of  anode  inlet  gas,  tda 

60 
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Fig.  6.  Capillary  effect  on  the  net  water  drag  a,  as  a  function  of  the  cathode  humidi¬ 
fication  temperature.  A:  calculation  without  capillary  effect;  a,  ♦,  and  ■:  calculation 
with  capillary  effect  (tceH  =  80°C,  tda  =  60°C,  £/=1.43,  £fl  =  2,  a:  pc,mpl  =  12 kPa,  ■: 
Pc,mpl  =  30  kPa,  ♦:  pc,mpl  =  48  kPa).  The  solid  vertical  bars  in  the  figure  represent  the 
experimental  data  distribution  ranges. 


Fig.  7.  Change  in  the  capillary  effect  on  the  net  drag  a,  as  a  function  of  the  cell  oper¬ 
ation  temperature  and  relative  humidity  (£/=1.43,  £a  =  2,  pc,mpl  =  301<Pa).  ■:  high 
humidity  condition  (RH  =  93-100%),  ♦:  low  humidity  condition  (RH  =  44-69%). 


exist  some  differences  in  the  porosity,  peak  pore  diameter  and 
surface  energy  of  these  carbon  papers  [12,13],  these  differences 
do  not  exceed  about  10%  of  the  values.  (2)  For  the  calculation  of 
capillary  pressure  characteristics,  it  is  necessary  to  know  both  the 
hydrophilic  and  the  total  PSD  of  a  GDL.  The  PSD  data  of  SGL24DC  by 
Kumbur  et  al.  is  the  only  available  data  which  gives  both  hydrophilic 
and  total  PSD. 

The  net  drag  was  calculated  for  varying  capillary  pressure  pc. 
The  capillary  pressure  of  the  MPL,  which  corresponds  to  sl2  =  0.056 
(calculation  process  shown  in  Section  2.4)  is  30kPa.  pc,mpl  is  var¬ 
ied  from  12  to  48  kPa  (a  range  of  ±18  kPa  was  selected  for  both  the 
upper  and  lower  limits  from  the  computed  value  of  30  kPa)  and 
the  results  are  shown  in  Fig.  6.  There  was  no  significant  difference 
in  the  computed  net  water  drag  for  all  pc,mpl  in  the  low  humid¬ 
ity  region  (tdc<65°C).  Flowever,  as  the  inlet  humidity  increases, 
the  net  drag  for  pc,mpl  =  30  and  48  kPa  decreased  abruptly,  which  is 
in  good  agreement  with  the  experimental  data  and  suggests  that 
the  hydrophobic  micropores  of  a  MPL  act  to  transfer  water  from 
the  cathode  to  anode  and  reduce  the  possibility  of  flooding  in  the 
cathode  GDL. 

Karan  et  al.  claimed  that  the  effect  of  the  MPL  does  not  signifi¬ 
cantly  increase  the  water  transport  from  the  cathode  to  the  anode, 
according  to  experimental  data  [10,11].  However,  the  effect  of  the 
MPL  is  prominent  under  high  humidity  conditions  and  the  effect  of 
the  MPL  on  the  net  drag  is  dependent  on  the  cell  operation  condi¬ 
tions.  For  example,  calculated  net  drags  for  a  GDL  without  MPL  and 
a  double-layered  GDL  (pc,mpl  =  30kPa)  are  shown  in  Fig.  7,  and  the 
cell  operation  conditions  are  the  same  as  those  listed  in  Table  3, 
except  the  cell  operation  temperature.  The  difference  in  the  net 
drag  is  larger  as  the  cell  operation  temperature  is  increased  and  as 
the  humidity  is  increased. 

4.  Discussion 


4.1.  Role  of  the  MPL  for  improvement  of  cell  performance 

The  MPL  has  another  role;  the  improvement  of  cell  performance. 
When  a  cell  is  operated  under  relatively  low  humidity  conditions,  it 
is  desirable  to  keep  the  MEA  wet  in  order  to  reduce  the  membrane 
resistance.  On  the  contrary,  when  a  cell  is  operated  in  the  saturated 
condition,  generated  water  should  be  quickly  removed  from  the 
reaction  sites. 

Cell  performance  was  computed  as  a  function  of  the  cathode 
humidification  temperature  and  the  results  are  shown  in  Fig.  8a. 
The  cell  operation  conditions  are  the  same  as  those  listed  in  Table  3, 
except  for  the  humidification  temperature  of  the  cathode  reactant 
gas.  For  all  regions,  the  cell  output  voltage  increased  with  the  capil¬ 
lary  pressure  of  the  MPL.  The  cell  performance  was  improved  with 
both  increases  in  the  cathode  humidification  temperature  and  the 
capillary  pressure,  and  the  results  suggest  an  optimum  capillary 
pressure  for  the  MPL  in  the  range  between  30  and  48  kPa. 

The  increase  of  the  cell  output  voltage  is  mainly  due  to  the  con¬ 
ductivity  reflected  by  the  water  content  X,  in  the  membrane.  The 
differences  in  cell  output  voltage  and  ohmic  loss  in  the  membrane 
with  and  without  the  MPL  ( pc,mpl  =  30  kPa)  at  several  cathode  inlet 
humidities  are  shown  in  Fig.  8b.  The  increase  in  cell  output  volt¬ 
age  is  almost  equal  to  the  decrease  of  ohmic  loss  of  the  membrane, 
except  under  low  humidification,  where  the  cathode  overpotential 
increases  significantly  due  to  the  local  current  density  distribution 
along  the  flow  channel. 


Fig.  8.  (a)  Effect  of  the  MPL  on  cell  performance  as  a  function  of  the  cathode 
humidification  temperature  (tce//  =  70 °C,  tda  =  60 °C,  £/=1.43,  £a  =  2,  ♦:  Pc,mpl  =  0,  ■  : 
Pc,mpl  =  12  kPa,  •:  pc,mpl  =  30kPa,  a:  pc,mpl  =  48  kPa).  (b)  The  differences  in  cell  out¬ 
put  voltage  and  ohmic  loss  in  the  membrane  with  and  without  the  MPL  at  several 
cathode  inlet  humidities  (tceu  =  70  °C,  tda  =  60  °C,  £/=  1.43,  =  2,  pc,mpl  =  30  kPa). 


4.2.  Relationship  of  MPL  and  the  catalyst  layer  structure 

When  the  pressure  of  water  transport  from  the  MPL  is  higher 
than  the  capillary  pressure  of  the  hydrophobic  pores  in  the  catalyst 
layer,  the  water  generated  in  the  cell  reaction  may  cause  flooding  in 


1854 


E.  Nishiyama,  T.  Murahashi  /  Journal  of  Power  Sources  196  (201 1 )  1847-1854 


the  catalyst  layer.  Therefore,  the  optimum  capillary  pressure,  i.e., 
the  optimum  PSD  for  the  MPL,  also  depends  on  the  PSD  and  its 
apparent  water  contact  angle  to  the  catalyst  layer. 

Uchida  et  al.  reported  that  the  PSD  for  the  catalyst  layer  consists 
of  primary  pores  (<0.04  p,m)  and  secondary  pores  (0.04-1  p,m),  and 
the  latter  corresponds  to  the  gas  paths  [22,23].  They  also  reported 
that  the  specific  pore  volume  distribution  of  the  best  performance 
cell  had  an  average  pore  radius  of  approximately  0.1  p,m.  The  capil¬ 
lary  pressure  of  the  hydrophobic  pores  with  a  characteristic  radius 
of  0.1  pan  is  approximately  425  kPa  for  a  contact  angle  of  0C  =  110°. 
When  the  capillary  pressure  of  the  MPL  (ranged  from  12  to  48  kPa) 
is  less  than  the  capillary  pressure  of  the  catalyst  layer,  then  flooding 
does  not  occur  in  the  catalyst  layer.  However,  more  information  is 
required  with  regard  to  the  contact  angle  and  pore  size  distribution 
in  the  catalyst  layer  to  obtain  a  certain  conclusion. 

5.  Conclusion 

Water  transport  through  the  gas  diffusion  media  of  a  PEMFC  was 
investigated  with  a  focus  on  the  role  of  the  MPL  on  the  cathode  GDL. 
The  capillary  pressures  of  the  MPL  and  GDL  were  estimated  using 
the  PSD  model.  PSD  functions  are  derived,  of  which  the  parame¬ 
ters  are  determined  by  fitting  to  the  measured  total  PSD  data.  The 
capillary  pressure  calculated  for  a  GDL  with  MPL  was  in  good  agree¬ 
ment  with  the  experimental  results  at  low  liquid  saturation.  For 
a  GDL  without  MPL,  the  calculated  capillary  pressure  was  also  in 
good  agreement  with  the  measured  capillary  pressure  under  water 
intrusion. 

Based  on  the  relationships  between  capillary  pressure  and  liq¬ 
uid  saturation,  mass  balance  equations  along  the  cathode  flow 
channel  were  solved  to  obtain  water  transport  rates  and  cell  volt¬ 
ages  for  various  feed  gas  humidification.  The  capillary  pressure  of 
the  MPL  relates  to  the  activity  of  water  at  the  cathode  catalyst 
layer/membrane  interface  and  has  the  effect  of  decreasing  the  net 
water  drag  from  the  anode  to  the  cathode.  The  net  drag  calculated 
from  the  PSD  model  with  consideration  of  the  capillary  pressure 
corresponded  well  with  the  measured  values. 

The  effect  of  the  MPL  capillary  pressure  on  the  cell  output 
voltage  was  significant  for  cell  operation  under  high  humidity 
conditions,  which  suggests  that  the  wetness  of  the  MEA  has  an  opti¬ 
mum  pore  size  distribution  and  the  contact  angle  of  the  catalyst 
layer  to  prevent  flooding. 
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